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We perform a numerical study of the preheating mechanism of particle production in models of 
quintessential inflation and compare it with the usual gravitational production mechanism. We find 
that even for a very small coupling between the inflaton field and a massless scalar field, g > 1CP 9 , 
preheating dominates over gravitational particle production. Reheating temperatures in the range 
1CT 2 GeV < T rh < 10 15 GeV can be easily obtained. 
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I. INTRODUCTION 

Inflationary models of the early universe have become 
a paradigm that can explain the isotropy of the cosmic 
microwave background, the origin of the small adiabatic 
perturbations that seed galaxy formation and the flatness 
of space-time. There is also evidence that the universe 
has recently entered another stage of mild inflation, re- 
sponsible for its accelerated expansion. It would be nat- 
ural to investigate models where the two periods of infla- 
tion are related. In particular, models of quintessential 
inflation explore the possibility that the same scalar field 
is responsible for both periods of inflation [1,2]. 

In quintessential inflation models, the inflaton poten- 
tial has no minimum in which the field could oscillate to 
produce the matter that would reheat the universe, as in 
usual chaotic inflation models. Therefore, it is assumed 
that all matter and energy in the universe will be created 
by gravitational production associated to changes in the 
geometry of the universe. However, it is known that this 
process is not very efficient [3]. In this letter we assess 
the importance of another mechanism of particle produc- 
tion in the early universe, namely preheating [4], in which 
particles are produced due to the variation of the clas- 
sical inflaton field. We find that the introduction of a 
non-zero coupling between the inflaton and matter fields 
leads to particle production that can easily overcome the 
gravitational one. 
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where A = 1 x 10~ 14 is required from structure formation 
[6] and M = 5 x 10 8 GeV in order to produce the present 
quintessential energy density [1]. 

Defining adimensional parameters <p = *&/Mpi, q = 
M/Mpi and r = VAMp;t, the classical equation of mo- 
tion for the inflaton field becomes independent of the A 
parameter: 
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where the primes are the derivatives with respect to r 
and we have used the Hubble parameter, 
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II. THE MODEL 

For definiteness, we study the Peebles and Vilenkin 
model [1], with a potential for the inflaton field $ of the 
form 



In order to study preheating, we couple the classical in- 
flaton field to a massless quantum scalar field \ through 
the lagrangian: 

c = \d^x - \iRx 2 - \g 2 <s> 2 x\ (4) 
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where R is the scalar curvature, £ is the gravitational cou- 
pling and g is a quartic coupling constant between <j> and 
X- Ford [3] calculated the energy density of x particles 
produced gravitationally in the limit |£ — 1/6 1 <C 1 due 
to the transition of the universe from the de Sitter space- 
time at the end of inflation. However, this production is 
not efficient since the ratio of matter and inflaton field 
energy densities is typically p m /p<f> ~ 10~ 18 [3]. After 
thermalization of the products it is possible to define a 
reheating temperature when radiation starts to dominate 
and in that case it is of the order of T rh <~ 10 4 GeV. 

When preheating is operative, the picture changes 
drastically. Taking into account the dilution of gravi- 
tational particle production, at the end of inflation, until 
the beginning of the particle production by the preheat- 
ing mechanism, the universe is still dominated by the 
inflaton density energy and this is our starting point. 

The equation of motion for each comoving k mode of 
the x ne ld can be written in terms of a new convenient 
variable X k = a 3 ^ 2 Xk as: 



III. PARTICLE PRODUCTION 



In order to compute the total amount of particles pro- 
duced during preheating, we have to solve equation (5). 
To obtain 4>(t) we choose the initial conditions as being 
4>{t = 0) = -1 and 4>'{t = 0) = 0, since in A$ 4 model 
inflation ends when $ ~ —Mp\. Then we can evaluate 
the evolution of the scale factor a(r) given by: 



a(r) = ao exp 



(It'H(t') 



(6) 



where H{r) is given in equation (3) and ao = 1. We 
show in figure 2 the evolution of the scale factor during 
the period relevant for particle production. 



(r) 



- A \X k = 0. (5) 



The physical momentum is given by k/a and we have 

verified that A = | (j^j + § (ir) i s neglegible during 
the production of \ particles. 

This is an equation of motion for an harmonic oscil- 
lator with variable frequency, lo\ = (k 2 / (Mp t a 2 (r)) + 
g 2 (f) 2 (t)) / \. When the adiabaticity condition is broken, 
that is, ut' k > lu 2 , particles in the k mode are produced. 
Figure 1 illustrates the time interval where adiabaticity is 
violated, leading to particle creation. This time interval 
corresponds to the passage of the inflaton field through 
its origin. 
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FIG. 2. Evolution of the scale factor as a function of r. 

Finally, substituting (6) and the solution of equations 
(1 — 2) into (5) we obtain X k (r). 

We define an adiabatic invariant n k that can be in- 
terpreted as the comoving number density of particles 
produced with momentum k: 



+ \Xk 



fe - o 



(7) 



FIG. 1. u>l (solid line), u' k (dotted line) and 10 x 
(dot-dashed line) as a function of t for j = 1CP 5 and 



k 



5vAMp;. Notice that adiabaticity is violated around 



t ~ 7, when the inflaton field goes through zero. 



In figure 3 we show the comoving number density of 
particles created during preheating obtained from the 
solution of equation (5) with vacuum initial conditions 
for different values of comoving momentum k, for g = 
5 x 10~ 5 . One can see that particle production occurs at 
t — 7, and quickly stabilizes. 



2 



M 



2 

1.75 
1.5 
1.25 
1 

0.75 
0.5 
0.25 



; 


jj|f*»— ^—^ 


^ k=100 ?i°- 5 M p i - 


V„ k=200 ?i°- 5 M pl * 









10 
x 



15 



20 



FIG. 3. Comoving number density of \ particles produced 
during preheating as a function of r for comoving momentum 
k = 0, 100 and 200 \f\M P i, for S = 5 x 10" 5 . 

In order to illustrate the fc— dependence of the den- 
sity of particles, we plot in figure 4 the density spectrum 
produced during preheating for different values of A and 
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FIG. 4. Comoving number density of % particles produced 
during preheating as a function of comoving momentum k in 
units of \f\M p . A = 1 x 10~ 14 and g = 5x 10~ 5 (solid curve); 
A = 1 x 10~ 14 and g = lx 10~ 4 (dashed curve); A = 4x 10~ 14 
and g = 5 x 10~ 5 (dot-dashed curve). 

The spectrum is well described by the expression: 

1 V\ k 2 



n k = exp 



52 g \M 2 pl 



(8) 



This result is to be compared to the approximation for 
particle production in the first oscillation of a pure (p 4 
potential [5]: 



n 
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= exp 
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XMpi k 2 /a 2 (T p ) 
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(9) 



where t p is the time of particle production and |<1> P | is the 
velocity of the inflaton field at the bottom of the poten- 
tial. In our case we have |$ p | = 3.6 x 10~ 2 \/\M PI and 
a(r p ) = 68, which results in: 
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exp 
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(10) 



an excellent approximation indeed. 

The total energy density produced at t p is given by: 
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f d 3 k 
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n k ~ 3 x lQ- e Xg 2 M Pl (11) 



Since the energy density in the inflaton field around 
the time of particle production is 



^ ~ 6 x 10" 4 AM^ , 



(12) 



the energy density that is transferred from the inflaton 
field to the \ ne ld & t the moment of production is given 
by: 



^~5x 10"V 



(13) 



This means that there is an upper limit in energy that can 
be extracted from <j> field, at the moment of production, 
~ 5 x 10~ 3 . However, once the x-particies are pro- 
duced they rapidly become non-relativistic, since their 
effective masses are given by m x = g<& (see equation(5)). 
In such a case, the above comparison becomes 



^ ~ 14 ff 5 / 2 , 



(14) 



for $ ~ 10 _1 Mp; which is a good approximation for early 
times as shown in figure 1. 

In order to prove the initial statement that the diluted 
energy density of particles produced gravitationally, p m 
[3] is in fact much smaller than the inflaton energy den- 
sity, at <p ~ 0, when the preheating starts, we take the 
ratio between them 1 : 



Pm((f> = Q) p m ((j) = ~l)/a~ 



P4,{<j> = 0) P4> (<j> = -l)/a- 6 
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-22 



(15) 



On the other hand, we can compare p m with our \ pro- 
duction (equation (11)): 



M0 = Q) 

Pm(<P = 0) 



10 23 ff 5/2 



(16) 



1 Although is not proportional to a 6 through all the in- 
terval, it corresponds to a maximum dilution. 
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and conclude that the latter dominates over the former 
for g > 6 x 1(T 10 . 

Finally, let us estimate the reheating temperature. We 
will make the assumption that the \ particles decay fast 
enough to avoid the effects of backreaction over the <j> 
field. In reference [5] it is shown that backreaction is 
important only for large times. So, without loss of gen- 
erality, we will consider the transference of x energy to 
relativistic particles just after the \ production. 

As the inflaton field energy density scales as a~ 6 since 
it is dominated by kinetic energy whereas the relativistic 
particles energy density, pu, scales as a~ 4 , the latter will 
start to dominate the universe at a time td given by: 



a 2 (T p ) 



Ug 5 / 2 



a 2 {r d ) 

At this time, the pa will be diluted to the value: 



(17) 



p R = 3 x 1(T 5 A 3/ V /2 ^pz $ 



a 4 (T p ) 



~ 6 x 10- 4 A 3 /y 5 / 2 M 4 ; , 

and the reheating temperature will be: 

T rh ~ (p R )^ T rh ~ KT 1 (X 3 ^g 16/2 ) 1/A M Pl . 

(18) 

For our limits on g, 6 x 10 10 < g < 1, the reheating 
temperature corresponds respectively to: 



10- 2 GeV < T rh < 10 15 GeV, 



(19) 



which is a confortable range. Further constraints in T r h 
imply in new bounds on the coupling constant between 
the inflaton and the x fields. 



IV. CONCLUSIONS 

Quintessential inflation models postulate that the same 
field is responsible for inflation in the early universe and 
the accelerated expansion of the universe today. In this 
type of models, the potential has no minimum and the 



traditional reheating mechanism of the universe can not 
be operative. Gravitational production of particles is 
usually adopted to reheat the universe. 

We have demonstrated in this letter that introducing 
a small coupling g > 10~ 9 between the inflaton field and 
a massless scalar particle x leads to particle creation in 
the preheating mechanism that dominates over the usual 
gravitational production in quintessential inflation mod- 
els. 

Reheating temperatures in the range 1CP 2 GeV < 
T rh < 10 15 GeV can be easily obtained for 6 x 1(T 10 < 

Our result is stronger than the one obtained in the ap- 
proximation used in [5] because we took into account the 
dilution of the produced particles between the epochs 
of gravitational creation and preheating. It is also not 
fully complete since the gravitational production mecha- 
nism could be affected when the x ~ 4> coupling is intro- 
duced in the lagrangian. However, the effect would be to 
inhibit the gravitational production, making our results 
even stronger. 
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